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ABSTRACT 

Recent measurements of the if— band luminosity function now provide us with strong, 
rehable constraints on the fraction of baryons which have cooled. Globally, this fraction 
is only about 5%, and there is no strong evidence that it is significantly higher in 
clusters. Without an effective sub-grid feedback prescription, the cooled gas fraction 
in any numerical simulation exceeds these observational constraints, and increases 
with increasing resolution. This compromises any discussion of galaxy and cluster 
properties based on results of simulations which include cooling but do not implement 
an effective feedback mechanism. 
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1 INTRODUCTION 

Gas cooling in the expanding Universe is an intrinsi- 
cally unstable process because cooling acts to increase 
the density of the gas, which in turn increases the cool- 
ing ra te. The consequence of this is that as soon as 



the matter in the Universe is contained within haloes. 



The best numerical models (Jenkins et al. 2000) in- 



the gai; within a system is able to cool at all it tends 
to do so catastrophically, only regulated by the speed 
at which the gas can respond to the new configura- 
tion. Systems which collapse at low redshift have a low 
mean density and, thus, long cooling times which gen- 
erally exceed their dynamical times. Systems collapsing 
at higher redshifts have a higher mean density and, be- 
cause the cooling rate depends strongly upon the gas 
density, they cool much more rapidly. In a Universe 
dominated by cold dark matter (CDM), the hierarchi- 
cal growth of structure consequently results in almost 
all of the baryons cooling by th e present day, unless 
additional p hysics is considered ( Cole 1991; White & 
Frcnk 199l|; [Blanchard et al. 1992 j; 



Unfortunately, the naive theoretical prediction that 
all the baryonic material in the Universe should have 
cooled into galaxies and formed stars conflicts with the 
observation of X-ray emission from galaxy clusters. This 
emission demonstrates that large amounts of hot gas 
persist in the Universe. This material can arise from two 
possible sources; either it was never contained within a 
collapsed halo in which cooling was efhcient, or energy 
injection (for instance feedback of energy due to super- 
novae) has either reheated or prevented the cooling of 
t he ga s. The commonly used theory of Press & Schecter 
(1974) (see also Bower 1991) implicitly assumes that all 



dicate that the fraction of material contained within 
haloes is certainly high, and only leave room for a small 
uncoUapsed f raction. Ther e fore, we appeal to fee d back 
mechanisms (Larson "197^; IWhite fc Rees 197^ ; |CoI^ 



1991; White & Frenk 1991) to reduce the amount of 



cold gas within collapsed haloes. 

In this paper we review the observational con- 
straints on the density of cold and hot baryons, both 
globally and in rich clusters (§^). These data establish 
the global cooled gas fraction at /c,giobai~ 0.073/i, and 
there is no convincing evidence that it is much higher 
in clusters. We then examine cosmological simulation 
results in §|[ and show that they can easily cool far 
too much gas to be consistent with these data. Finally, 
in §^ we discuss the necessity of combining simulations 
with analytic, sub-grid physics models, and the merits 
of various initial attempts to do this. 

Throughout this paper, we parametrise the Hubble 
constant as Ho = IQOh km s^^ Mpc^^. We will need to 
consider two cosmologies, A CDM (fio — 0.3, A — 0.7) 
and standard CDM (fio = 1, A = 0). Matter densities 
are given relative to the critical density, and mass-to- 
light ratios {M/L) are in solar units. 



2 OBSERVATIONAL CONSTRAINTS 

2.1 The Global Fraction of Cold Baryons 

The mass density of cooled baryons in the form of stars, 
^^stars, can be obtained from the observed luminosity 
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density (e.g. Blanchard, Valls-Gabaud & Mamon 1992) 
as long as 1) we do not miss a significant contribution 
from low surface brightness galaxies; and 2) the mass- 
to-light ratio is well understood. Fortunately, surface 
brightness is tightly correlated with total luminosity, 
and there does not appear to be a significant contribu- 
tion to the baryon density from low surface brightness 
galaxies (Driver 1999; Cole et aL2000). Although there 
is an uncomfortably large range in the Schechter param- 
eters of the local optical luminosity function (due to the 
large d egeneracy between (tf and L*), Fukugita, Hogan 



CO-based compilation of Young & Scoville (1991). To- 



& Pee bles (1998, hereafter FHP) demonstrate that the 



total integrated _B— band luminosity is constrained to 
within about 15%. However, the use of the optical lu- 
minosity function requires large and uncertain correc- 
tions for the M/Lb ratio, which is sensitive to stellar 
populations and, thus, must be determined separately 
for different types of galaxy. 

A more robust estimate is obtained from the 
i^— band luminosity function, since M/Lk is less sensi- 
tive to star formation history, varying by only a factor 
of about two due to stellar population ag e differ ences 
( iBellfcde Jong 20001 ). Recently, Cole et al. ( |2000bD have 
computed the mass function using over 17000 galaxies 
with redshifts from the 2dF galaxy redshift survey and 
magnitudes from the 2-micron all sky survey (2MASS). 
From the galaxy colours and stellar population synthe- 
sis models they determine M/Lk for each galaxy (inde- 
pendent of h). Their {AC DM) result is most sensitive 
to the assum ed initial mass function (IMF): for a Ken- 
nicutt (I1983D I MF, fJ^tars = 0.0014/1"!, ^^iHe for the 
Salpetei~( |r955| ) IMF, f^^tars = 0.0026ft." ^ 

Th e average stellar M/Lk from which Cole et al. 
(2000b) derive rJstars is 0.73 and 1.32 for the Kcnni- 



cutt and Salpeter IMFs, respectively (independent of 
h). A substantially higher M/Lk would require a large 
contribution from low mass stars and brown dwarfs. Al- 
though these objects are likely to be a significant com- 
ponent of the mass budget, there is no strong evidence 
that they ar c more abundant than expected from stan- 
dard I MFs (|Fuchs et al. 1998] pizis et al 20001 ; [Lucas 
& Roc re 200D ). In particular, Reid et a/. ( 199!j ) claim 
that brown dwarfs contribute no more than 15% of the 
mass of the Galactic disk. Earlier microlensing results 
which suggested that there might be a large population 
of dark objects in the Galaxy halo are no longer com- 
pelling, following a reanalysis of the MACHO project 
data (Alcock et al. 2000). 

Dynamical measu r ements of early type galaxies 
( Ivan der Marel 199l| ; |BeU fc de Jong 2000| ) favour 
M/Lk <j l.l/i; this is an upper limit to the global 
average because later type galaxies have lower M/Lk- 
We therefore prefer the estimate of rJstars — 0.0014/1"^ 
based on the Kcnnicutt (1983) IMF with a mean 
M/Lk = 0.73. 

Cooled gas is also present in the form of neutral and 
molecular gas, but these make only small contributions. 
The density of neutral gas is riatomic ~ 0.000188/1"^, 
as recently calculated by the HI Parkes Al l Sky Sur- 
vey (HIPASS) team ( |Kilborn et al. 19991 ). For the 



gether, the neutral and molecular gas only represent 
~ 10% of the stellar mass, which is the fraction we will 

adopt, so that rJcold = l-lfigtars- 

The total baryon content of the universe is unknown 
observationally, because baryons in the warm plasma 
phase believed t o occupy normal gala xy haloes are very 
difficult to see ( Benson et al. 2000| ). However, if nu- 
cleosynthesis calculations are correct, the baryon den- 
sity rib can be determined from deuterium abun dances. 
The best current estimate, Vlh = 0.019/i~^ ( Buries 
fc Tytler 1998^ [Buries et al. 200C| ), therefore imphes 
/c. global = f^coid/r^b = 0.073/i. Thcrc arc indications, 
from recent analysis of the combined BOOMERANG 
and MAXIMA data, that Of, is much larger than this, 
Vtb = 0.039/1-2 pg^ffg aL2000), which would reduce 
/c,giobai by a factor of two. 

In summary, only about 5% of the available baryons 
in the universe have cooled (for h k, 0.7 and a Kennicutt 
IMF). This result is slightly lower than the minimum of 
the range 0.062 < /c,giobai < 0.167 found by FHP from 
optical data, but much more certain as it does not rely 
on knowing the relative abundances of different galaxies 
(disks, spheroids and irregulars) with very different (and 
uncertain) M/Lb ratios. As we will show in §||, this 
low value is in stark contrast to the results obtained in 
numerical simulations which do not employ an adequate 
feedback model. 



2.2 The Abundance of Hot Baryons 

In the previous section we claimed that most of the 
baryons in the Universe are in a warm or hot phase, 
which is difficult to observe because it is not hot enough 



to be observable in X-ray radiation (Blanchard et 



1992; Dave et al. 2000; Benson et al. 200C). The ex 



istence of such a warm component is compatible with 
constraints on the anisotropy of the microwave back- 
ground, a long as the gas temperat ure satisfies the con - 
straint T < 4/i x 10^ K (at z < 1) ( [Wright et al. 1994] ). 

In clusters, however, the same gas is much hotter, 
and is direct ly observable. As originally shown by White 



et al. (1993), the total baryon content of rich clusters, 
including this plasma, is fully consistent with fif, from 
element abundance determinations, if Vto ~ 0.3. FHP 
compile data on the baryon contributions in clusters 
due to stars, cold gas and plasma, and demonstrate 
that 9.b/VLo = 0.112 ± 0.05 (/i = 0.7). In this case, 
ilfc is a sum of all observed baryons in clusters. Hence, 
for 0.0 = 0.3 ± 0.1, we have that Ob = 0.034 ± 0.019 
{h = 0.7) which is consistent with the baryon fraction 
dedu ced from deuterium abundances (Buries fc Tytler 
1998| ). Thus, the fact that fistars ^ in the Universe 



amount of gas in molecular form, we adopt the rela- 
tion PH2/PHI = 0.81 used by FHP, gleaned from the 



is almost certainly a reflection of the fact that the warm 
gas has not been directly observed. 



2.3 The Dependence on Halo Mass 

Although we have demonstrated that /c,giobai~ 0.073/i, 
it is possible that the cooled fraction in some environ- 
ments, rich clusters for example, could be very differ- 
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ent. In particular, Bryan ( 2000D claims that the efh- 
ciency of galaxy formation is dependent on halo mass, 
with more gas cooling in small haloes, relative to large 
ones. Most of the cluster data in the literature has been 
obtained at optical wavelengths, so the correction for 
M/L variations is not as robust as for the global value 
c om puted from the if-band field luminosity function in 
§ pn| . We will use the optimal values from FHP, who 
compile data from various sources, including dynami- 
cal measurements and population synthesis model es- 
timates. The M/Lb depends strongly on morpholog- 
ical type, and FHP find M/Lb = e.Sla.o ^r E/SO 
galaxies; M/Lb = 1-5 ± 0.4 for spiral galaxies; and 
M/Lb = 1.1 ± 0.25 for irregular galaxies. For the mor- 
phological composition of a "typical" cluster, they de- 
termine M/Lb = 4.5 ± 1, and we will use this number. 
Assuming a pure elliptical population would increase 
stellar mass estimates by ^ 45%. 

For low mass clusters and groups, the X-ray data 
only extend out to a small fraction of the virial radius. 
This can lead to a strong bias, since th e gas distribu- 



-^^coid and Mhot are distance dependent measurements, 
with Mcoid oc h^^ and Mhot oc h^^-^ . Since, usually, 
M\^ot ^ Afcoid, /cool hP'^ when calculated from Equa- 
tion 0. 

Most of the /c, cluster values in the bottom panel 
of Figure |^ are higher than the global fraction /c, global 
computed in § |2.1| , shown as the horizontal line. There is 
only weak evidence for a trend, in the sense that /c, global 
is lower in higher temperature systems. Howeve r, this 
trend is much weaker than suggested by Bryan ( 2000[ ) 
and, we caution that the kT < 5 keV groups may still 
be biased toward high fractions because of the limited 
extent of the X-ray data; none of these groups have 
X-ray emission detected beyond half the virial radius. 
Also, both the random and systematic uncertainties in 
observationally determining Mhot, especially for the low 
temperature groups, are large. This is emphasized in 
the top panel of Figure |l|, in which the same data are 
plotted, but /c, cluster is computed as 



tion is less concentrated than the stars ( Roussel et al 
200(]()7]We will consider the cluste r and group data re- 



fc. 



Mcold 



:lustc 



total 



(2) 



cently compiled by Roussell et al. ( 2000 ) specifically for 
the purpose of addressing this issue. In that paper, care 
was taken to ensure that stellar and gas masses are com- 
puted to the same radius, and only those clusters for 
which the X-ray data extend to at least 25% of the 
virial radius are considered here. Nonetheless, we warn 
that the data for groups with kT < 5 keV must be 
treated with caution, as their gas masses still require a 
substantial extrapolation to the virial radius. In addi- 
tion to the standard analysis of the X-ray data, Rous- 
sel et al. also consider hot gas mass estimates based on 
simulation-motivated scaling laws, which do not require 
an assumption of isothermal, /3— model profiles. We use 
the measurements made based on the scaling laws of 



Bryan & Norman (1998). In addition, we consider the 
15 rich clusters with re 



where Mtotai is the total dynamical mass of th e clus - 
ter, and we take fit from Buries & Tytler ( |l998[) . 
This does not depend on the difficult measurement of 
Mhot, but does require an estimate of fio, which is 
0.3 in our ACDM cosmology. Note that, in this case, 
/c.ciusteroc floh. Not Only has the apparent trend in 
/c. cluster with temperature now virtually disappeared, 
but the scatter has been reduced considerably. In this 
case the /c. cluster values are in much closer agreement 
with fc. ,giobai; tli6 remaining difference is not very com- 
pelling, since a slightly lower Qo ~ 0.2 would result in 
concordance between the cluster and global frac tions. 

These results differ from those of Bryan ( ^000 ), 
who claims that /c, cluster increases from cluster to group 
scales. However, this trend could well be the result of 
a bias due to the limited radial extent of the X-ray 



iable velocity dispersions from obser vations, i?^, in the group samples ( Roussel _et 



the CNOC sample (lOarlberg et al. 1996; KJarlbcr g et a l 



r g et a l 



i997|)7Fmr A-ray data taken trom Lewis et a, 

The stellar masses are computed from the inte- 
grated stellar luminosity function, assuming M/Lb = 
4.5, and corrected for undetected galaxies by extrap- 
olating the luminosity function with an assumed faint 
end slope of a « —1.2. Finally, we assume that the mass 
of cold baryons is 10% larger than the stellar mass, to 
accou nt fo r the presence of neutral and molecular gas 
(see § 2.1). 

We consider two approaches to compute /c, cluster 
from A/cold in the clusters. In the ffi'st case, we compute 

/c, cluster — Mcold/ (Mcold + Mhot), (1) 

where Mhot is determined from the observable X-ray 
emission from the hot plasma, extrapolated to the virial 
radius. This will be valid if the extrapolation is accu- 
rate (less likely for lower mass clusters), and will be 
an overestimate of /c, cluster if there is an unobserved, 
warm component which does not contribute to the X- 
ray emission (as is the case in galactic mass haloes) . We 
show /c,ciiister for the two samples in the bottom panel 
of Figure ^ for a ACDM cosmology with h ~ 0.7. Both 



200C). For example, from the Mulc haey et 
sample considered by Bryan ( pOOO ), the groups wit 



the smallest have among the highest /c, duster values. 
On the other hand, the two groups for which > 0.5 
Mpc have the lowest /c, cluster- This emphasizes the im- 
portance of ensuring the X-ray data extend out to at 
least a sizable fraction of the virial radius, as we have 
done when considering the Roussel et al. compilation. 
We have shown that there is no discernable trend of 
/c, cluster with temperature in this sample, particularly 
if the baryon fraction is assumed from nucleosynthesis 
arguments, rather than relying on the measurement of 
Mhot. 



3 NUMERICAL SIMULATIONS 

Hydrodynamical numerical simulations which attempt 
to follow the effects of gas cooling have proved very diffi- 
cult to perform across the full range of the gravitational 
mass hierarchy. Allowing gas to dissipate its energy via 
radiative cooling generates enormous density contrasts 
and a huge range of interesting spatial scales within the 
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Figure 1. The fraction of baryons in tlie cold phase in clus- 
ters, as a function of X-ray temperature. The data are from Carl- 
berg et al. (1996, circles) and Roussel et al. (2000, squares), la 
error bars are only available for the Carlberg et al. sample. The 
cluster data are all normalized to a stellar M/Lb = 4.5; increas- 
ing this number increases /c, cluster proportionally. The horizontal 
line shows the best estimate of the global cold baryon fraction, 
/c, global I based on the X— band luminosity function results of 
Cole et al. (2000); note that global h. Bottom panel: /c, duster 
is computed from Equation hi, where the total baryon mass is 
the sum of observed galaxies and hot gas. The data points scale 
approximately like /<^iustor « h°-^ . Top pane/./c^dustor is com- 
puted from Equation H, which does not depend on a measurement 
of Miiot, but is dependent on f2o. In this case, duster « Q.oh. 



model, even before the intrinsically sub-resolution pro- 
cesses of star formation and the feedback of energy due 
to supernovae and stellar winds are considered. 

The earliest attempt at s uch a simulation was made 
by Thomas & Couchman ( 1992 ), though the model 
which included radiative cooling was only briefly dis- 
cussed. This model had a very high mass resolution 
threshold, allowing only the largest galaxies to form, 
and only a small fraction of the gas to cool. Th is was 
quickly followed by the work of Katz & White (1993) 
who produced a model that was several years ahead of 
its time. They demonstrated that the formation of a 
Virgo-like cluster, including the effects of cooling, could 
be followed by a computer simulation. However, the 
model had some problems: in particular, as the authors 
make clear, much more gas cools than is observed. These 
results are shown as the large, filled circles in Figure 
0. The lower point shows the actual fraction of cold 
baryons in the simulation; the higher point to which it 
is connected shows this fraction corrected by extrapo- 
lating the simulated luminosity function to include the 
additional cooled gas within unresolved objects, as de- 
scribed in Katz & White. The simulation of Katz & 
White was stopped at z — 0.13, so the redshift zero 
fraction will be even higher. 



Figure 2. The solid points are the cooled baryon fractions in 
several standard CDM simulations, as labelled. The lower point 
refers to the measured cooled fraction, and the upper point to 
which it is joined is corrected to include the contribution of unre- 
solved galaxies. All of these simulations use /t = 0.5, and cannot 
be rescaled for other values of h. The Dave et al. (2000) work 
considers both treecode and Eularian grid-based simulations in 
several cosmologies; we represent the average cooled fraction of 
these simulations as a lower limit, because it does not include 
cooled gas in unresolved galaxies. The data are shown as open 

fmbols, and are the same as those in the bottom panel of Figure 
but for h = 0.5 instead of h = 0.7. 



Recently there has been rei iewed interest in this 
problem. Suginohara & Ostriker ( 199^ ) also produced a 
cluster model, and noted that additional physical pro- 
cesses (such as star formation and feedback) would be 
required to prevent the cooling of an excessively large 
fraction of the gas; unfortunately, the fraction of cooled 
baryons was not published in their paper so we can- 
not consider their results in Figure B. Their work was 



foll owed by that of Pearce et al. ( |1999| ; |2000| ), Lewis et 
(2000) and Dave et aZ. (2000). Lewis et al. used a 
more advanced code to im prove upon the earlier sim- 
ulation of Katz & White ( 1993| ), and they also found 
that around 40% of the gas within the cluster cools. 
This is shown in Figure as filled triangles; the up- 
per triangle includes a correction for galaxies below the 
resolution limit, by extrapolating the resultant lumi- 
nosity function assuming a faint end slope a = —0.96. 
Dave et al. used both a parallel treecode and a Eulerian 
grid based code in several cosmologies and, although the 
cluster-by-cluster fractions are not quoted, find a uni- 
versal cooling fraction of around 30% to 40% for their 
high resolution simulations, in both cases. This is shown 
as a lower limit in Figure since it does not include a 
correction for galaxies formed in haloes below the reso- 
lution limit. 

Simulations like those discussed above which cool 
too much gas cannot be expected to model galaxy 
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formation in a physically correct way. An alternative, 
heuristic approach can be taken to force the simulations 
to satisfy the cons traint on /c. global, by limiting the res- 
olutio n (see § pTl] ). This is the method used by Pearce 
et al. ( ^000 ), who deliberately set their mass resolution 
to cool ~15% of the gas within their simulation volume; 
the cooled fractions for the twenty most massive of these 
clusters are shown in Figure as the filled squares (both 
before and after correction ior sub-resolution galaxies). 
By construction, the cooled gas fractions of these sim- 
ulated haloes are closer to the observed values, though 
they are still probably too high. The obvious problem 
with this simulation is that it does not allow the forma- 
tion of small galaxies. 

3.1 Resolution 

It has long been theorised that progessively increas- 
ing the resolution within a simulation leads to a con- 
tinual increase in the fraction of material which cools 



(Cole 1991; White fc Frcnk 1991; Suginohara fc Ostrikei 



199S), a situation commonly referred to as the cooling 
catastrophe. We demonstrate this in Figure ^, by pre- 
senting the results of standard cold dark matter simula- 
tions, at different resolutions (the mass of the smallest 
object that can effectively cool). This figure shows how 
the fraction of gas particles in the cold phase increases 
as resolution is increased. At the low resolution end, the 
relation is too steep, due to the artifici al he ating prob- 
lems described by Steinmetz & White (1997). However, 
the trends discussed below are qual itative y r etained, 
even when this problem is overcome (Kay 2000). Apart 
from this effect, the shape of the curve will be approx- 
imately independent of the simulations. However, the 
scale along the abscissa will vary; that is, the curves 
may be translated laterally, for example, by changing 
r^b. More subtly, in a model with a positive cosmological 
constant, structure forms earlier and so more gas cools 
at the same threshold mass. This causes the curves plot- 
ted on Figure |^ to move to the right by about a factor of 
two. Changes in metallicity can affect the result in the 
same way, though even primorid al abundances can not 
prevent the cooling catastrophe (Dave et al. 2000). In 
addition, differences in the numerical implementation 
can also shift the curve by factors of order two. 

Figure ^ makes it clear why some groups find cooled 
gas fractions around 40% whilst others find much lower 
fractions. Values around 40% will be achieved for a wide 
range of simulation resolutions, as this is the value ob- 
tained when the resolution mass is less than the char- 
acteristic mass of the halo mass function. More gas will 
cool as the resolution is increased but a large increase 
is required to see a significant change, since the slope of 
the mass function is shallow in this regime. Any re- 
quired cold gas fraction below 40% can be obtained 
by choosing the resolution appropriately, although care 



0.6 - 



,0.4 - 



0.2 - 



Shocked 



Galaxy 



Uncollapsed 

A , , , , A 



10'° 10" 
min. baryon mass of galaxy/Mg 

Figure 3. Fraction of baryons in three phases (representing 
shocked gas, cooled gas [galaxies] and unshocked, diffuse baryons) 
against the baryonic mass of the smallest resolvable galaxy for 
standard CDM simulations with radiative cooling. The simula- 
tions use a box-size of lQh~^ Mpc and f2b = 0.06. The two points 
with the largest mass threshold are too low, due to artificial heat- 
ing effects. Further details are contained in Kay (2000). 



"cooling crisis" ; it is a crisis because its resolution will 
result in a significant improvement in the utility of sim- 
ulation results. 



4 DISCUSSION 

We have demonstrated that it is easy to exceed the ob- 
served cooled gas fraction with numerical simulations 
of sufficient resolution. The failure to match this con- 
straint reflects the well-known need for some method 
of reheating the gas in ord er to preven t excessive star 
form a tion at high redshift ( Larson 1974| ; White & Rees 



must be taken to ac count for artificial heating (Stein- 



metz & White 1997). However, in order to get a good 
match to the observed gas fraction within a standard 
cold dark matter simulation the resolution limit must 
be set to ~ 10^^ Mq or higher, which is clearly unphys- 

ically high. This is the situation we referred to as the of this gas is cold and dense, and so any excess energy 



1978| ; [White fc Frcnk 1991| ). Until an effective feedback 
model is implemented in simulations, it will not be pos- 
sible to obtain physically illuminating results. 

Present attempts to incorporate feedback within 
cosmological simulations all suffer from the same draw- 
back, that the physical processes responsible for the re- 
distribution of energy occur on scales well below those 
that can presently be resolved. Hence, an accurate treat- 
ment of the propagation of energy from the local ISM 
to super-galactic scales is currently an intractable prob- 
lem. However, several groups have attempted to include 
feedback by injecting the energy into regions at the reso- 
lution limit of the simulation. For example, an approach 
that has been taken in smoothed particle hydrodynam- 
ics (SPH) simulations is simply to add available energy 
from newly- formed "star" par ticles to the thermal en- 
ergy of the surrounding gas (Katz et al. 1996 ; Lewis 
' 20001) . However, the majority 



et al. 2000; Dave et 
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is rapidly re-radiated. It is plausible that this problem 
is due to the failure of the SPH simulations to resolve 
a multiphase medium, in which the feedback would be 
able to create hot bubbles of diffuse gas, with signif- 
icantly longer cooling times than the cold dense gas 
from which they were spawned. Several groups have at- 
tempted to "correct" for this by preventing reheated 
gas from radiatively cooling over a finite timescale (e.g. 
Gerritsen 1997; Springel 2000; Thacker & Couchman 
2000). Another method is to supply energy in kinetic 
form (e.g. Navarro & White 1993; Kay 2000), which is 
efficient at limiting the fraction of cooled gas. However, 
the effectiveness of this implementation may be arti- 
ficial since the SPH simulations fail to resolve shocks 
t hat would efficien tly thermalize the reheated material 
( |Katz et al. 1996[ ). 

An alter nativ e approach to feedback is that of Cen 
& Ostriker ( |1999D , who use a grid-based cosmological 
code, and implement an alytic rules, similar to those 
of semi-analytic models (Bomcrville fc Primack 1999|; 



Cole et al. 2000a; Kauffniann et al. 1999| ; Somerville 
et al. ^OOip to regulate the formation of galaxies within 
a single cell. Unfortunately, the resolution of grid-based 
methods such as this is still too coarse (~ 100/i~^ kpc 
in the highest resolution simulation of Cen & Ostriker) 
to provide useful results on galaxy scales. It is not clear 
that the analytic rules implemented to model the sub- 
grid physics in this code are robust to changes in res- 
olution, as changes in the cell size can result in large 
changes in the gas densities, for which the analytic 
formalism may not appropriately compensate. A more 
promising route may be the use of adaptive-mesh rc- 
finement codes (Bryan & Norman 1998; Bryan & Voit 



200l|)7]which will allow the implementation of an ana- 
lytic treatment of gas cooling within cells of sufhcient 
resolution. 

The goal now is to find a physically motivated feed- 
back model which is able to sufficiently reduce /c. It is 
therefore useful to consider, qualitatively, what is re- 
quired of such a model. The amount of gas which can 
cool in a given halo cannot be computed by consider- 
ing a halo structure at a single redshift and treating 
it as if it had existed unchanged over a Hubble time. 
For example, the cooling time in clusters today is much 
longer than the Hubble time, and this has been used 
to explain w hy not all of the gas in t h ese environments 
has cooled (Binney 1977; 3ilk 1977; Rees & Ostrikci 
1977|)~| However, this argument docs not hold in hi- 
erarchical models, where cluster progenitors at earlier 
epochs were low mass haloes in which the cooling time 
was short, and simulations of hierarchical cluster for- 



Mcooi, as the mass which can cool all the gas out to its 
virial radius, within a Hubble time at redshift z. The 
cooling rate, which depends on halo density and tem- 
perature, is computed using the models of Raymond et 
al. (1976), with one-third solar metallicity. This cooling 
mass, for the standard CDM model described above, is 
shown in Figure ^, and is given by Mcooi ~ 10^^ M©, 
almost completely independent of redshift for z < 10. 
This non-trivial result is a consequence of the fact that, 
when cooling is dominated by line emission, the increase 
in the cooling rate with increasing redshift is closely bal- 
anced by the decrease in the Hubble time. This forms 
an approximate upper mass limit to efficient cooling; 
haloes with masses smaller than this limit can cool all 
of their gas, out to the virial radius, while to the right 
of this line, only a rapidly decreasing fraction of the gas 
is able to cool. 

At low redshifts, most of the haloes which are col- 
lapsing are more massive than Mcooi; therefore, cooling 
plays a minor role today. By z k, \, however, the char- 
acteristic mass crosses the Afcooi threshold, and cool- 
ing becomes a dominant process, leading to the cool- 
ing catastrophe. In simulations, the resolution thresh- 
old imposes a mass limit, Mth, which prevents smaller 
mass haloes from cooling. Efficient coolingis therefore 
restricted to the shaded region in Figure ^. As Mth is 
lowered, more and more mass cools, and will quickly ex- 
ceed the tight observational constraints. If the threshold 
is high enough to prevent efficient cooling at high red- 
shift, no low mass galaxies are formed. 

The final fraction of cooled gas then depends on the 
trajectory a particular halo takes through this plane, as 
it grows. This growth will occur via both smooth accre- 
tion, and discrete jumps due to mergers. The resolution 
threshold imprints a scale on the accumulation of cooled 
baryons, as the final fraction /c in a given halo will de- 
pend on its trajectory only after it becomes sufficiently 
resolved. This could have an effect on the dependence 
of /c on halo mass in simulations (but see Kay 2000). 

It is clear that any analytic feedback model which 
behaves in a manner similar to a sharp, fixed mass 
threshold will fail to simultaneously match the observed 
/c, global and still produce a reasonable galaxy luminos- 
ity function. A physical feedback scheme, however, does 
not provide a hard mass threshold, but gradually alters 
the cooling efficiency of haloes as a function of mass 
and redshift; thus the shaded region in Figure I be- 
comes a "fuzzy" region without a well defined lower 
boundary. There are many models for s uch feedback , 
including supernovae or AGN reheatin g ( Larson 1974 ■ 
White fc Rees 19781; |Laccy fc Silk 199T|; |Cole 1991|; [Wii 



et al 19981; [Bower et al. 2000| p>/Iac Low &: Ferrara 
1999| ; [Efstathiou 2000| )_^ ionizing radiat i on from the ul- 



matioE demonstrate that a large fraction of the baryons 



cool bj the present day. To treat the problem analyti- 
cally, i t is necessary to consider the hierarchical growth 



of a halo through the plane of Figure 4. The smooth 
curved lines represent the evolution of structure, for a 



travi olet background (Efstathiou 1992; Weinberg et 
1997 ) , and an entropy threshold established at high rcd- 



standa pTCDM model (cts ^ 0-8, P ^ f^p/i ^ 0.5); the 
central line shows the median mass halo which virial- 
izes as a function of redshift, and the two others bracket 
the region in which 50% of the mass in the Universe col- 
lapses int o halo es. In a manner similar to that of Rees & 
Ostriker (1977), we can define an efficient cooling mass. 



shift ( Kaiser 1991|; Blanchard et al. 1992 ; Ponman et 



1999| ; balogh et al. 19991 ). A very successful implemen 



tation of a feedback scaling law is t hat used by semi- 



analytic models of galaxy for n iation ( ^ole et al. 2000a 



Somerville & Primack 1999; Kauffmann et 



1999 



Somerville et al. 2001), in which the mass of reheated 



gas is assumed to scale in proportion with the halo cir- 
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Figure 4. The region of efficient gas cooling in the Universe, as 
a function of halo mass and redshift. The smooth curved lines 
delineate the formation of structure in the Universe. The central 
line represents the median halo mass which virialises at the cor- 
responding redshift; 50% of the mass in the Universe collapses 
between the two lines which bracket this one. Mcool is the mass 
threshold below which haloes can cool out to their virial radius. 
A/th is the mass threshold within the numerical model, generally 
given by the resolution, below which cooling cannot occur. The 
shaded region therefore denotes the region within which cooling 
is important. 



cular velocity raised to some power a. Though this scal- 
ing is constructed to allow the models to match the ob- 
served luminosity function, it also succeeds in limiting 
/c to more reasonable values (though generally still high 
compared with the constraints described here). Further- 
more, there may be some theoretical motivation for this 
scaling relatio n from supernova-driven wind models (Ef- 
stathio u 200C ). The next step in this work is to combine 
these models with realistic halo merger trees to deter- 
mine, analytically, how they affect the cooling rate as a 
function of halo mass and redshift, for implementation 
into numerical simulations. 



5 CONCLUSIONS 

In this paper, we have presented an up-to-date review 
of the cooling catastrophe. Although the issues we have 
discussed have been known for a long time, recent ad- 
vances in observations and simulation techniques war- 
rant a reevaluation of the situation, which can be sum- 
marized in two key points: 

• The global fraction of baryons in the form of stars 
and cold gas is now well constrained observationally 
from the if— band luminosity function to be /c,giobai~ 
0.073/i. Although there is some uncertainty about how 



this may vary with halo mass, there is no convincing ev- 
idence for a significant difference from the global value. 

• The amount of gas which cools in simulations is 
directly related to the resolution, unless a working re- 
heating scheme is implemented. No current simulation 
is able to produce both (1) the correct /c, global and (2) a 
population of low mass galaxies, without implementing 
analytic models of the sub-grid physics on uncomfort- 
ably large scales. 

The use of simulations to aid our understanding 
of galaxy properties and of cluster scaling laws there- 
fore awaits the implementation of a sub-grid, physi- 
cal feedback scheme which allows the models to sat- 
isfy observational constraints on /c, global in a resolution- 
independent way. 
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